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Abstract 
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1 Introduction 



Lattice QCD simulations employing a discretisation of the Dirac operator 
based on the original proposal by Wilson [l] are currently being performed 
with light dynamical fermions [2]j9]. When decreasing the light quark mass 
at a fixed value of the lattice spacing, a subtle interplay between mass and 
discretisation effects can take place due to the explicit breaking of chiral sym- 
metry by the Wilson term. In simulations with light values of the quark mass 
it is, therefore, vital to understand and monitor the discretisation effects and 
to obtain a quantitative measure of their size. 

Close to the continuum limit, a useful way to determine the discretisation 
effects in the regime of light quark masses is provided by Wilson chiral pertur- 
bation theory (W^PT), an extension of the continuum chiral effective theory 



including additional terms proportional to powers of the lattice spacing 10 
Depending on the order of the expansion, additional low energy constants 
(LECs) appear, whose values are not known a priori: they depend on the 
lattice action and can only be computed from a simulation. 

Knowing the values of the LECs of W%PT is of particular interest, because 
W%PT predicts a non-trivial phase structure for Wilson type fermions in the 



lattice spacing and quark mass plane 10 -12 . Depending on the sign of a par- 



ticular combination of Wilson LECs - commonly denoted as c-i cx — (2Wg+Wg) 



- either the Aoki-scenario 11 or a first order, so called Sharpe-Singleton [10] 
scenario is realised. Numerical evidence for both scenarios has been observed 
in lattice QCD simulations and dedicated studies of the associated phase di- 



agrams have been performed by several groups 13 -24 



The Aoki scenario with a positive C2 was found to be realised in quenched 
simulations. In dynamical simulations the Sharpe-Singleton scenario with neg- 
ative C2 was observed when using Wilson twisted mass fermions at maximal 
twist 20 23 . This manifests itself in the fact that the neutral pion mass M n o is 
lighter than the charged one, M n ± , where the splitting of the squared masses is 
proportional to c^- In turn, a measurement of the pion mass-splitting in Wilson 
twisted mass lattice QCD provides a way to measure C2 and, hence, the LECs 
of the corresponding chiral effective theory. However, this way of computing 
C2 is challenging since for the neutral pion mass disconnected contributions 
need to be evaluated. 



There are also alternative ways to determine the LECs of W%PT. In Refs. 



28 



34 



they have been studied by matching the analytical predictions 27 



25- 



29- 



for the spectrum of the Wilson Dirac operator - with fixed index in a 
finite volume - to lattice dataJZ] Determinations of the Wilson LECs have 



For a recent review, we refer to 135 



2 



36 38 



also been carried out via the spectral density of the Hermitian Wilson-Dirac 
Lattice determinations of the pion scattering lengths have 
an approach that was extended to a partially 
. In a mixed action with Wilson-type sea fermions 



operator 
been used in Refs. 39l 40 



41 



quenched setup in Ref. 
and chirally invariant valence quarks, a mixed action chiral Lagrangian 42143 



can be constructed. The corresponding LECs - in particular W$ - have been 
recently determined in the case of overlap valence quarks on a Wilson twisted 
mass sea 



44 



In this paper we are going to determine the LECs using a method recently 
proposed in Ref. 45 . It relies on the measurement of pseudoscalar meson 



masses involving Wilson twisted mass fermions. In this approach, the LEC 
Wg is related to the mass-splitting between the charged pion mass M n ± and 
the "connected neutral pion" mass M^(o,c) . The latter is determined from the 
quark-connected correlation, which contributes to the neutral pion in twisted 
mass QCD. The "connected neutral pion" correlation function thus differs 
from the complete correlation function needed to determine the neutral pion 
by the absence of disconnected diagrams. This has the advantage that in 
numerical studies the "connected pion mass-splitting" should be accessible 
with good statistical precision. In addition, the mass-splitting, M^ (0 c) — M^ , 
between the connected and full neutral pion mass, provides an estimate of W§. 
We are also going to study the mass and lattice spacing dependence of these 
splittings arising at higher order in the W%PT expansion. 



The paper is structured as follows: in the next section we collect the W%PT 
expressions relating the pion mass-splittings to the Wilson LECs. In Section [3] 
we present the lattice actions used in our study. The determination of the 
LECs from two of those lattice setups is described in Section |4j Finally, a 
qualitative comparison of the values of these LECs from different choices of 
the lattice action is reported in Section [5j 



2 Wilson Chiral Perturbation Theory (W%PT) 



In this section, we briefly discuss Wilson chiral perturbation theory and intro- 
duce the expressions used in our study. For a recent review on the applications 



of ;\RT to lattice QCD, we refer to Ref. 46 



Our study will be based on the computation of pseudoscalar meson masses 
involving Wilson twisted mass fermions. We are therefore interested in a chiral 
Lagrangian involving a mass matrix M = + ififr^, where and fi^ are 
the renormalised untwisted and twisted quark masses, respectively. The masses 



m and /j,£ appear in the Wilson twisted mass action of eq. (12). 
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At leading order (LO) in the power counting, m ~ fie ~ a 2 AQ CD , and after 
a shift of the quark mass to remove a term of 0(a), the partially quenched 



chiral Lagrangian reads 10 43 



C x = L str (d^Zd^) - f- 8 ^ Str (il^E + E f M 



Str (E - E f 



-a 2 W^ [Str (E + E f )] 2 - a 2 W^ 
-aXStr(E 2 + [E t ] 2 ) , (1) 

where E parametrises the vacuum manifold and thus characterises the Nambu- 
Goldstone bosons arising from the spontaneous breaking of chiral symmetry. 
The traces over the flavour indices are denoted by Str and a = 2W a. In 
addition to the continuum LECs Bq and / (defined with the convention giving 
f n ~ 130 MeV), the chiral Lagrangian also includes Wo, W$, Wj and W%, which 
are Wilson LECs describing discretisation effects. 

In this work we are interested in the determination of the Wilson LECs by 
matching lattice QCD calculations of pseudoscalar meson masses to their 
PQW^PT expressions. As already mentioned, we consider Wilson twisted 
mass fermions at maximal twist. This is achieved in the chiral Lagrangian 
by setting mo = 0. At non vanishing values of the lattice spacing, the break- 



ing of flavour symmetry by the twisted mass term in eq. (12) implies that the 
mass of the charged pion M n ± differs from that of the neutral pion M T o by 
0(a 2 ) effects. A similar pattern holds for the mass of the "connected neutral 
. The PQW^PT expressions for these three meson masses 



(0,c) 



pion mass, 
at LO read f45l|47l449 



M 2 ± 
M 2 



t(0 ,c) = 2B ^ e - 8a 2 w' 8 , 
where w' k is related to the Wilson LEC W' k by 



Ml 



(2) 
(3) 
(4) 



im 2 w> 

w k — 



P 



(k = 6,8) 



(5) 



To determine the individual values of the LECs it is, therefore, possible to 
consider the following mass-splittings: 



M 2 ± 



Ml 



(0,c) 



,2 „../ 



■5a w 



8 • 



\ (M 2 (0 , c) - M 2 ) = 8 



a 2 w'p. 



(6) 
(7) 
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From eq. ([3]), it appears that the linear combination of LECs which controls 
the mass-splitting between charged and neutral pions is given by 

19W 2 

C2 = -^(2WZ + W2). (8) 
This can also be re-expressed as 



c 2 



-2(2w' 6 + w' 8 ) 



(9) 
(10) 



In the W%PT expressions presented above, two light mass-degenerate flavours 
were assumed to be present in the sea sector. When considering also other 
dynamical flavours, such as the strange and the charm quarks, the same ex- 
pressions hold when assuming that these heavier flavours sufficiently decouple 
from the light quark sector. In this case, the values of the Wilson LECs will 
have a further residual dependence on the heavier quark masses. 

Before closing this section, we mention that W%PT calculations at NLO have 



been carried out for the pion mass and decay constant with Nf = 2 50-52 
and Nf = 2 + 1 + 1 [531 flavours of twisted mass fermions. 



3 Lattice actions 



The complete lattice action can be written as 

S = Sf+Sg, (11) 

where Sf is the fermionic action and S g is the pure gauge action. As we shall 
see, in this work we will consider a few alternatives for both the fermionic and 
the gauge actions in order to explore the dependence of the Wilson LECs on 
the details of the lattice action. As discussed below, we will use for Sp a few 
variants of Wilson twisted mass fermions. 



3.1 Wilson Twisted Mass Fermions 



The Wilson twisted mass (Wtm) lattice action for the mass degenerate light 

(12) 



doublet (u, d) in the so-called twisted basis reads 54 , 55 
Si 



° 4 iXi(x) [D[U] + m ,i + i/J,tj 5 T 3 ) Xi(x)} , 
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where m ,z is the untwisted bare quark mass, fie is the bare twisted light quark 
mass, r 3 is the third Pauli matrix acting in flavour space and 



Dp] = \ [i, (v„ + v;) - av;vj , (13) 

is the massless Wilson-Dirac operator. V M and V* are the forward and back- 
ward gauge covariant difference operators, respectively. Twisted mass light 
fermions are said to be at maximal twist if the bare untwisted mass m 0) z is 
tuned to its critical value, m crit . The quark doublet xi — (Xu, Xd) m the twisted 
basis is related by a chiral rotation to the quark doublet in the physical basis 

^Phys = e ^ n5 r 3 ^ ^phys = - t e ^ n5T3 ^ 
where ui is the twist angle. 

The twisted mass parameter provides an infrared regulator avoiding the 
presence of accidental zero-modes in the Wilson-Dirac operator. An important 
property of Wtm fermions is that at maximal twist physical observables are 
0(a) improved |55|. In numerical simulations, maximal twist is achieved by 
tuning the value of the hopping parameter k = l/(2rao,z + 8) to its critical 
value K crit by tuning the PCAC quark mass m PC Ac to zero. The expected 
0(a 2 ) scaling of physical observables when performing the continuum limit 



extrapolation has been confirmed in the quenched approximation 56 59 and 



with Nf = 2 60 63 and Nf = 2 + 1 + 1 (7 64 65 dynamical quarks. 
A peculiar lattice artifact can appear in observables made out of Wtm quarks 



due to the breaking of isospin and parity by the twisted mass term in eq. ( 12 ). 
This effect, which is expected to vanish in physical quantities at a rate of 0(a 2 ) 
when approaching the continuum limit, has been observed to be numerically 
small in most of the observables which have been analysed [7,60-64 



An exception to this observed small isospin breaking effects is found in the 
case of the neutral pseudoscalar mass. Indeed, isospin breaking induces a 
mass-splitting between charged and neutral pion masses. Dedicated numer- 
ical studies indicate that while in the charged pion mass only very mild cutoff 
effects are present, the neutral pion mass is instead affected by significant 
C(a 2 ) effects [7j[60j[63j[65}[66]. An analysis based on the Symanzik expansion 
indicates that isospin breaking affects only a limited set of observables in a 
sizeable way, namely the neutral pion mass and kinematically related quanti- 
ties 67,68 . This analysis is complementary to that based on W%PT where, as 
previously mentioned, the mass-splitting between charged and neutral pions 
is parametrised by the combination of LECs appearing in C2 defined in eq. (pi). 



The determination of the neutral pion mass M n o involves both connected and 
disconnected contributions. The computation of quark-disconnected diagrams 
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is challenging and requires the employment of specific techniques in order to 
achieve a statistically significant determination of M n o 69 70 . 



N{ — 2 + 1 + 1 Wtm fermions 

In addition to a doublet of mass-degenerate light quarks (u,d), a heavier dou- 
blet with strange and charm quarks - (s,c) - can be incorporated in lat- 



tice QCD studied with Wilson twisted mass fermions 71 72 . Also in this 
Nf — 2 + 1 + 1 setup, a tuning to maximal twist by imposing mpcAc = 0, 
allows to achieve the automatic 0(a) improvement of physical observables. We 
refer to Ref. \7\ for a complete description of the lattice action for the heavier 
quark doublet and for further details on the lattice setup. 



Mixed action with Osterw alder- Seiler valence quarks 



Osterwalder-Seiler (OS) valence quarks 73 can be viewed as the building 
blocks of Wilson twisted mass fermions at maximal twist. The OS action for 
an individual quark flavour x / reads 

Sf = a 4 E {*/(*) [D[U] + ™ ciit + i» f <y 5 r f ] X f{x)} , (15) 

X 

where rj (here |r/| = 1) is the Wilson parameter and m cr it the critical mass. By 
combining two flavours of OS quarks with opposite signs of r/, e.g. r 2 = — r 1; 
the action of a doublet of maximally Wtm fermions of mass \ij — \i\ — can 
be recovered. The benefits of the OS action are that 0(a) improved physical 
observables |71| can be obtained by using the same estimates of m crit as in 
the Wtm case and, thus, avoiding further tuning effort. OS and Wtm fermions 
coincide with Wilson fermions in the massless limit and consequently share the 
same renormalisation factors. This simplifies the matching of sea and valence 
quark masses in the context of a mixed-action with Wtm sea and OS valence 
quarks. 

The pseudoscalar correlation function obtained when considering only the con- 
nected contribution to the neutral pion correlation function (i.e. when ignor- 
ing disconnected diagrams), is precisely the pion correlator with OS fermions. 
The mass-splittings in eqs.(|6])-([7]) can hence be interpreted as involving sea 
and valence quarks in a mixed action setup. 

In this work, we aim at determining the Wilson LECs in a lattice theory with 
Nf = 2 + 1 + 1 Wtm fermions. A particular effort will be devoted to addressing 
the main systematic effects present in these determinations. We furthermore 
aim at exploring the qualitative change on the values of these LECs when 
varying the details of the lattice action. Below, we briefly summarise the al- 
ternative lattice setups used in this work. We will consider variants of the 
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action differing by the presence of smearing of the gauge links in the covariant 
derivative, the inclusion of the Sheikholeslami-Wohlert term or by a change in 
the number of dynamical flavours. 



Stout Smearing. 



A smearing procedure can be applied to the gauge links entering in the covari- 



ant derivatives in eq. (13). The stout smearing 74 procedure is analytic in 



the un-smeared link variables and hence well suited for simulations with the 
HMC algorithm. The smearing can be iterated several times, with the price 
of extending the coupling of fermions to the gauge links over a larger region. 



Sheikholeslami- Wohlert term. 



In our comparison of the values of the mass-splittings in eqs. Q-0 from differ- 
ent lattice setups, we will also consider results available in the literature from 
quenched lattice simulations with Wtm quarks including the Sheikholeslami- 
Wohlert term 



75 



3.2 Gauge action 



The lattice gauge actions considered in this work have a generic form which 
includes a plaquette term U^ u and rectangular (1x2) Wilson loops U^ 2 U , 



^ = fE( & o E {1-BbTt(U™ )} + !>! £{1 



ReTr(^fJ} 



(16) 



fj,,i/=i 

l<fl<V 



with (3 = 6/<7q the bare inverse coupling and the normalisation condition 
b = 1 — 8&i. We will consider the case of the Wilson plaquette [I] action 
(&i = 0), the tree-level Symanzik improved 76 77 action {pi = —1/12) and 
the Iwasaki 78-80 action (b\ = —0.331). In Wtm simulations, the strength 



of the phase transition has been found [21 22 to depend on the value of the 
parameter bi in eq. (16). 
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4 Numerical Studies 



4-1 iVf = 2 + 1 + 1 Wtm fermions with Iwasaki gauge action 



The purpose of this study is to determine the mass-splittings in eqs. @ and 
([7]), which are directly related to the Wilson LECs W' & and Wq, respectively 
The lattice action is composed of the Iwasaki gauge action and Nf — 2 + 1 + 1 
flavours of Wilson twisted mass fermions. 

The simulations [7] were performed at three values of the lattice gauge cou- 
pling, (3 = 1.90, 1.95 and (3 = 2.10, corresponding to values of the lattice 
spacing a m 0.09 fm, 0.08 fm and 0.06 fm, respectively. The charged pion mass 
M n ± approximately ranges from 230 MeV to 510 MeV. Simulated volumes cor- 
respond to values of M n ±L larger than 3.3. Physical spatial volumes range from 
(1.9 fm) 3 to (2.8 fm) 3 . 

The values of the pseudoscalar meson masses for the Nf = 2 + 1 + 1 ensem- 
bles are collected in Table [l] The mass-splittings in eqs. g-Q, which are 
directly proportional to W$ and W' 6 respectively, are illustrated in Fig. [TJ We 
observe that the conditions W% < and W$ > are fulfilled by the lattice 



data, in agreement with the bounds derived in Refs. 29 , 34 , 45 . In order to 
quote the values of the Wilson LECs W$ 8 , the systematic effects from quark- 
mass dependence, residual lattice artifacts and finite volume effects have to 
be addressed. 

In the "large cut-off effects" power counting and at LO in the W%PT chiral 
Lagrangian, the mass-splittings in eqs. (|6|-((7]) are expected to be independent 
of the lattice spacing and the light-quark mass. The possible presence of such 
effects might thus signal effects entering at higher orders in the W%PT chiral 
expansion. Since the NLO expressions for these mass-splittings is currently 
not available in the literature, we rely in our systematic error analysis on a 
separate study of (i) the continuum-limit of the mass-splittings at a reference 
mass and (ii) the comparison of a constant and a linear chiral extrapolation 
in Ml±. 

Starting with point (i), we show in Fig. [I] the mass splittings eq. (|6| and ^ 
relevant for W$ and W$, respectively, as a function of M 2 ±. We first observe 
that data points with similar values of (M n ±r ) 2 , but coming from different 
lattice spacings, tend to be compatible with each other, in particular when 
considering the larger lattice sizes, represented by the filled symbols. The 
lattice spacing dependence of the two aforementioned mass splittings at a 



reference mass (M v ±r ) 2 pa 0.55 is illustrated in Fig. 2(a) Note that only the 
largest lattice sizes L are considered in this figure. Although the lattice size 
slightly varies when changing (3, the lattice data fulfils L > 2.5 fm and M W ±L > 
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Table 1 

Determination of the charged M w ± , neutral connected M V ( , C ) and neutral M n o pseu- 
doscalar meson masses from simulations with Nf = 2 + 1 + 1 flavours of Wilson 
twisted mass fermions at maximal twist and the Iwasaki gauge action. The value 
of the Sommer scale tq is determined in the chiral limit. The ensemble names con- 
taining an "s" or a "c" refer to a change of the strange or the charm quark mass, 
respectively. The values for M n ± agree within errors with the once published already 
inRef. (7|. 

4 and, therefore, we do not expect a large effect from a small mismatch in the 



physical volume. Fig. 2(a) suggests that the residual lattice spacing effects are 
small. We remind that these lattice artifacts appear beyond the leading 0(a 2 ) 
effects. As an estimate of these effects, we include in our systematic error 
analysis the difference between the values of the mass-splittings from the two 
finer lattice spacings. 

Concerning the light-quark mass dependence in point (ii), we can expect that 
the mass terms appearing at NLO can contain a linear term in M^ ± but 
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(b) 



Fig. 1. Determination of the mass-splittings (a) (M^ ± — M 2 {0c) )/a 2 in eq. (j6j) and 
(b) (M 2 (0jC ) — M 2 )/2a 2 in eq. (7h as a function of M 2 ± . These quantities are in 
units of the chirally extrapolatecTSommer scale ro- The mass-splittings are directly 
related to the Wilson LECs W$ and Wq. The lattice setup with the Iwasaki gauge 
action and A^f = 2 + 1 + 1 flavours of Wilson twisted mass fermions is considered. 
Filled and empty symbols signal a change in the lattice size. The label "(s,c)" in 
the legend indicates the effect of changing the values of the strange and charm 
quark masses in the sea. The mass-splittings illustrated in these figures assess the 
size of 0(a 2 ) discretisation effects. Hence, any lattice spacing dependence in these 
quantities points to residual higher order discretisation effects. The results of a chiral 
extrapolation from a constant and a linear fit in M 2 ± are shown. The deviation 
between the two extrapolated values is included in the systematic error analysis. 

also a term of the form M 2 ± log(M 2 ±). Indeed, such terms are present in the 
W%PT expression of the mass splitting between the charged and neutral pion 
masses at NLO 



51 



Since the precise form of the these logarithmic terms is 
yet unknown for the mass-splittings considered here, we limit ourselves to a 
linear chiral extrapolation in M 2 ± . Note that our data is not precise enough to 
disentangle possible logarithmic contributions. We take as our central values 
the linearly extrapolated mass-splittings and use the difference with respect 
to the constant fit as an estimate of the systematic error from the light-quark 
mass dependence. 



Finite volume effects are taken into account by adding to the systematic error 
the difference between the values of the mass splittings from two ensembles - 
A40.24 and A40.32 - differing only by a change of lattice size from L w 2.1 fm 
to 2.8 fm. This is expected to be a conservative choice since these ensembles 
correspond to a rather small light-quark mass - and therefore finite size effects 
can be non- negligible. Also, these ensembles were obtained at the coarsest 
lattice spacing, where possible finite size effects (FSE) from the neutral pion 
mass should be larger. 
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Fig. 2. (a) Lattice spacing dependence of the mass-splittings |M 2 ± — M 2 (0c) |/a 2 
- empty symbols - and (M 2 (0c) — M 2 )/2a 2 - filled symbols - at a reference 
mass (M^-iro) 2 ~ 0.55. For a better visibility, the absolute value of the mass 
splitting is considered in the case of the empty symbols, (b) Pion mass-splitting, 
(M 2 — M 2 ± ) /a 2 , normalised according to eq. Q in order to relate it to the com- 
bination of Wilson LECs C2. The chiral extrapolation using a constant and a linear 
fit in M 2 ± is also shown. The lattice setup with the Iwasaki gauge action and 
Nf = 2 + 1 + 1 flavours of Wilson twisted mass fermions is considered. 



As already mentioned, the determination of the Wilson LECs W' & 8 from lattice 
data with Nf — 2 + 1 + 1 flavours assumes that the strange and the charm 
sea-quarks decouple sufficiently from the light-quark dynamics. The residual 
heavy quark mass dependence present in W'§ 8 can be studied by varying the 
strange and charm quark masses in the neighbourhood of their physical values. 
This effect is illustrated in Fig. [l]by the points labelled by "(s,c)" in the legend. 
We use the difference between the mass-splittings from ensembles A80.24 and 
A80.24s - which only differ by a change of the strange and charm quark masses 
- to estimate this systematic effect. We expect that this is a conservative choice 
because (a) the ensemble A80.24s has a strange quark mass which is very close 
to the physical point, (b) the change in the strange quark mass is largest for 
A80.24 and A80.24s and (c) the effect of strange sea-quarks should be larger 
than that of charm quarks. 



After combining the previously discussed systematic uncertainties in quadra- 
ture, we obtain the following values for the mass-splittings for the case of a 
lattice setup with Nf = 2 + 1 + 1 Wtm fermions and the Iwasaki gauge action, 
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Table 2 

Determination of the Wilson LECs Wg 8 (u> 6 §) and C2 from a lattice setup with 
iVf = 2 + 1 + 1 Wtm fermions and the Iwasaki gauge action. For the values quoted 
in physical units, the input = 0.45(2) fm has been used. The values in the last 
column derive from eq. ([5]), where the input rof = 0.276(12) from [7] has been used. 
For the case of Wq 8 the systematic error analysis described in the text has been 
incorporated in the overall uncertainty indicated by the label "syst." in the table. 
For C2, in the last block, we quote separately the results of a constant and a linear 
chiral extrapolation in M 2 ± - as illustrated in Fig 2(b) - and the quoted errors are 
purely statistical. 

( Ml ± - M 2 , 0c) \ , , N 

' n * (0 =1 r £ = -23.0 ±0.7 ±3.0, (17) 



a 2 
2a 2 



ri = + 13.8 ±0.6 ±5.6, (It 



— -T 1U.U _L_ W.U _L_ U.U , 



where the first error is statistical and the second systematic. The correspond- 
ing values of the Wilson LECs are collected in Tab. [2} As already anticipated, 
the results for w' 6 8 are precise enough to identify a definite sign for these LECs. 



The combination of LECs C2, can be determined directly from the mass- 
splitting M 2 - Ml ± as indicated in eq. d9J). The measurements of this mass- 



splitting are illustrated in Fig. 2(b) where the results of a chiral extrapolation 
by using a constant and a linear fit in M 2 ± are also shown. For the case of 
C2, we provide in Tab. [2] the results from both these chiral extrapolations and 
quote in the individual numbers only the statistical error. The values arising 
from these extrapolations are both compatible with a negative sign of C2. 



The W%PT expressions at NLO relevant for C2 have been derived in Ref 51 
In addition to the Wilson LECs appearing at LO and to the usual Gasser- 
Leutwyler LECs, other parameters also appear at NLO. A complete determi- 
nation of these parameters is beyond the scope of this study. We postpone such 
an analysis to a future dedicated study of the W^PT description of both the 
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pion mass and decay constant. First results for the Gasser-Leutwyler LECs, 
from fits based on continuum xPT, have been presented in Refs. (7, 81 . 



The connected neutral pion mass, M w (o,c) , is an important ingredient in order 
to isolate the individual values of the LECs W$ 8 . As already pointed out, the 
connected neutral pion can be interpreted as the pion of a mixed action with 
OS fermions. Such a mixed action has been used to determine observables in 
the Kaon sector 



65,82,83 



We note that extensions of the analytical expres- 
sions to SU(3) W%PT is currently not available in the literature. Contrary to 
the pion case, the absence of disconnected diagrams in correlations functions 
in the Kaon sector could possibly allow to consider quantities from which the 
Wilson LECs can be determined with good accuracy. 



4-2 Nf = 2 Wtm fermions with tlSym gauge action 



In this section, we again determine the LECs W§ 8 , but this time using Nf — 2 
flavours of Wilson twisted mass fermions and the tree-level Symanzik improved 
gauge action [3, 63 70 . We already anticipate that a smaller set of ensembles 



and of measurements of the relevant pion masses are available in this case, 
in comparison to the Nf = 2 + 1 + 1 case discussed previously. Therefore, 
the resulting determinations and comparisons might suffer from insufficient 
control of systematic effects. However, we think that already a qualitative 
comparison can provide useful information to parametrise the size of cutoff 
effects from different lattice setups. 



The simulations considered in this work 63,85 were performed at three values 
of the lattice gauge coupling (3 = 3.90, 4.05 and (3 = 4.20, corresponding to 
values of the lattice spacing a ~ 0.08 fm, 0.07fm and 0.05 fm, respectively. 
The charged pion mass M n ± approximately ranges from 310 MeV to 460 MeV. 
Physical spatial volumes range from (1.3 fm) 3 to (2.6 fm) 3 and ensembles which 
differ only by the lattice size have been considered in order to address the size 
of finite volume effects in the determination of the Wilson LECs. 



The values of the pseudoscalar meson masses [63] f or the Nf = 2 ensembles are 
collected in Table[3] The mass-splittings in eqs. ([6])-([7]) are illustrated in Fig. [3j 
In order to explore the systematic effects present in these determinations, we 
follow a similar path to that described for the case of Nf = 2 + 1 + 1 ensembles. 
The availability of ensembles differing only by the physical volume allows to 
address the size of FSE in the mass-splittings. This is illustrated in Fig. |4(a)[ 
At P — 3.90, a set of three ensembles with L/a = 16, 24 and 32 could be used 
for the case of (M^± — M^ (0c) )/a 2 . With the current statistical uncertainties, 
no clear signs of FSE can be observed in the data. Furthermore, data from two 
different lattice spacings - f3 = 3.90 and 4.05 - agree within errors, indicating 
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Table 3 

Determination of the charged M w ± , neutral connected M n ( , c ) and neutral M w o pseu- 
doscalar meson masses from simulations with Nf = 2 flavours of Wilson twisted mass 



fermions at maximal twist and the tree-level Symanzik improved gauge action 63 
The value of the Sommer scale ro is determined in the chiral limit |84|. 
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Fig. 3. Determination of the mass-splittings (a) (M 2 ± 

(M 2 {0c) — M 2 ) / 2a 2 in eq. (|7|) as a function of M 2 ± . These quantities are in units of 



Q as a function of M 2 
the chirally extrapolated Sommer scale ro- The mass-splittings are directly related 
to the Wilson LECs W% and Wq. The lattice setup with the tree- level Symanzik 
improved (tlSym) gauge action and iVf = 2 flavours of Wilson twisted mass fermions 
is considered. Filled and empty symbols signal a change in the lattice size. 
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Fig. 4. (a) Finite volume effects on the mass-splittings \M^± — M 2 (0c) |/a 2 - empty 
symbols - and (M 2 (0c) — M 2 )/2a 2 ~ filled symbols. Note that the absolute value 
of the mass-splitting is used in the case of empty symbols, (b) Pion mass-split- 
ting, (M 2 — M 2 ± )/a 2 , normalised according to eq. (9j) in order to relate it to the 
combination of Wilson LECs C2- The chiral extrapolation using a constant and a 
linear fit in M 2 ± is also shown. The lattice setup with the tree-level Symanzik im- 
proved (tlSym) gauge action and Nf = 2 flavours of Wilson twisted mass fermions 
is considered. 

the absence of large residual lattice artifacts in these mass-splittings. However, 
the lack of sufficient data does not allow to address the mass dependence of 
the mass-splitting (M 2 (0iC) — M 2 )/a 2 in a satisfactory way. In analogy to the 
Nf — 2 + 1 + 1 case, we include the deviation between a constant and a linear 
extrapolation in M 2 ± in the estimate of the systematic uncertainties. The 
central value is taken from the result of the linear fit. 

For the case of a lattice setup with Nf = 2 Wtm fermions at maximal twist and 
the tlSym gauge action, we obtain the following values for the mass-splittings 



'K^-KA 4 

v 2a 2 



- 20.1 ±2.3 ± 1.7, 
+ 8.4 ±3.3 ±5.5, 



(19) 
(20) 



where the first error is statistical and the second systematic. The correspond- 
ing values of the Wilson LECs are collected in Tab.[4j The LEC W$ has recently 
been determined from a mixed action involving the same Nf = 2 lattice action 
in the sea sector as that described here, but with Neuberger overlap valence 

-0.0064(24), differs 



quark 44 . The value quoted in Ref. 44 



from the estimate in Tab. [4] at the 2-sigma level. However, we stress once more 
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Table 4 

Determination of the Wilson LECs Wg 8 (w' 6 8 ) from a lattice setup with Nf = 2 
Wtm fermions and the tlSym gauge action. For the values quoted in physical units, 
the input tq = 0.45(2) fm has been used. The values in the last column derive from 
eq. ([5]), where the input rof = 0.275(6) from [63] has been used. 

that our present Nf = 2 estimate does not include a complete assessment of 
the systematic errors. 



The determination of c 2 from Nf = 2 ensembles is illustrated in Fig. 4(b) More 
data would be needed to isolate the residual mass-dependence present in c 2 . 
For this reason we opt for quoting separately the results of a constant and a 
linear chiral extrapolation in M^± , 



C 2 To [est.] = 

c 2 Tq [lin. ] = 



-3.1 + 0.4. 
-0.5 + 1.5. 



(21) 
(22) 



where the errors are statistical only. These values are consistent with those 
arising from the measurements of pseudoscalar meson masses in Ref. [63] and 
are also very similar to the result obtained in Ref. 



50 for K 



-4c2 from 

twisted mass finite volume effects. All the lattice measurements favour a neg- 
ative sign of c 2 . However, as already mentioned, more data would be needed 
to properly address the residual mass dependence. 



We refer to Refs. 63 , 84 for more details about the description of the pion 



mass and decay constant by means of x?T expressions including discretisation 
effects. 



5 Discussion 



In this section we collect a few comments concerning the extraction of the 
Wilson LECs. The LECs W§ and c 2 depend on the neutral pion mass M n o. 
The statistical error in M w o is dominated by the contribution of disconnected 
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diagrams]^ We observe that the precise form of the light-quark mass depen- 



dence of the mass-splittings related to Wg and c 2 - see e.g. Figs. |l(b) and 



2(b) - cannot be addressed within the present uncertainties. As previously 
discussed this mass dependence can arise at NLO in the W%PT expansion. 
We can, however, not exclude at this stage that this higher-order effects are 
negligible. This issue is particularly relevant for the case of c 2 oc — (2w' 6 + w' s ), 
where a partial cancellation of the effect of w' 6 and w' 8 is present. The presence 
of higher order effects in the determination of c 2 has also been discussed in 
Ref. |40|. We recall that the sign of c 2 controls the appearance of an Aoki 
(c 2 > 0) or of a Sharpe-Singleton (c 2 < 0) scenario for the phase structure of 
Wilson fermions. 

From the previously discussed determinations, a comparison of the values of 
the Wilson LECs from the lattice actions (i) Nf — 2 + 1 + 1 flavours of Wtm 
fermions at maximal twist and Iwasaki gauge action and (ii) Nf = 2 flavours 
of Wtm quarks and tlSym gauge action, suggests that Wq 8 and c 2 do not 
vary significantly in between these two setups. Let us extend this observation 
by performing a comparison of the pion mass-splittings as determined from 
different lattice actions. We stress that a complete assessment of the overall 
uncertainty is not available for most of these measurements and therefore the 
comparison remains at the qualitative level. 

Lattice simulations with Nf = 2 + 1 + 1 Wtm fermions and the Iwasaki gauge 
action, but including in addition one iteration of stout smearing - labelled 
1-stout - have been reported in Ref. A qualitative comparison of the effect 
of the stout smearing on the size of the mass-splittings is shown in Fig. [5] Note 
that a single ensemble is used in the estimate of the mass-splittings for the case 
of stout smearing. The used smearing seems to help in reducing the magnitude 
of the splitting (M 2 ± — M 2 (0iC) ) / a 2 while - given the current uncertainties -it 
does not introduce a significant change in (M 2 (0jO ) — M 2 )/a 2 . 

Fig. [5^a) also includes the estimates of the mass-splitting (M 2 ± — M 2 (Q c) ) / a 2 
as determined from quenched ensembles either with or without the presence of 
the Sheikholeslami-Wohlert term. Maximally twisted-mass fermions and the 
plaquette gauge action are used in both cases. The values of the mass-splitting 
are derived from studies available in the literature. For the case without the 



Sheikholeslami-Wohlert term, results are taken from Ref. 66 . We use the 
information from different lattice spacings and quark masses to estimate the 
size of the systematic effects in (M 2 ± — M 2 (0c) )/a 2 . For the case in which the 



Sheikholeslami-Wohlert term is included - labelled csw ~~ we follow Ref. 86 



2 It is interesting to note that Table [T] indicates that the relative error on the neutral 
pion mass is roughly independent of the light-quark mass and that it decreases when 
increasing the volume. In practice, in the current simulation conditions, this implies 
that the measurements of M w o are statistically more precise for the ensembles with 
lighter quark masses. 
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Fig. 5. Comparison of the values of the mass-splittings (a) (M 2 ± — M 2 (0c) )/a 2 and 
(b) (M 2 ( 0c) — M 2 )/2a 2 from different lattice setups (see text for details). 

where the non-perturbative determination of csw was used. 

For cases other than those involving the Sheikholeslami-Wohlert term or stout- 



smearing, Fig. 5(a) suggests that, with the current precision, the value of the 
mass-splitting (M 2 ± — M 2 (0c) )/a 2 does not significantly depend on a simul- 
taneous change of the number of flavours N{ and of the gauge action. Note 
however that we cannot exclude that a change only in Nf or only in the pa- 
rameter hi of the gauge action, leads to a different conclusion. It would be 
very desirable, if further actions are investigated and the precision could be 
increased. 



86 is that 



One important observation arising from the measurements in Ref. 
the introduction of the Sheikholeslami-Wohlert term in quenched studies sig- 
nificantly reduces the size of 0(a 2 ) effects by lowering the value of the mass- 
splitting (M 2 ± — M 2 (0jC ))/a 2 . It would thus be interesting to study whether 
this result still holds for simulations with dynamical fermions and whether 
the value of (M 2 (0 c) — M 2 )/2a 2 is also reduced in that case. 



Conclusions 



We have presented the determination of the Wilson LECs Wq 8 and c 2 , para- 
metrising the size of 0(a 2 ) lattice artifacts in W%PT, from simulations with 
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a lattice action composed out of the Iwasaki gauge action and iVf = 2 + 1 + 1 
flavours of Wilson twisted mass fermions at maximal twist. The values of 
Wq 8 include a rather complete account of the systematic uncertainties. Our 



measurements satisfy the recently derived bounds 29 , 34 45 , W% < and 



Wq > 0. 

We have also explored the dependence of the mass-splittings (M 2 ± — M 2 (0jC) ) / a 2 
and (M 2 (o iC ) — M 2 )/2a 2 on the choice of the lattice action. From this quali- 
tative comparison, it is tempting to conjecture that a lattice action made of 
dynamical twisted mass fermions including the Sheikholeslami-Wohlert and 
smearing might lead to a reduction of the mass-splitting (M 2 ± — M 2 (0c) )/a 2 . 
Further studies are needed to clarify this point and to extend it to the case of 
(MVo - M 2 )/2a 2 . 

A partial cancellation of the contributions from W§ and W$ implies that the 
residual mass-dependence of c 2 is more sensitive to higher order effects in 
the W%PT expansion. While this potential reduction of c 2 would certainly be 
beneficial, the precise determination of its value cannot be achieved with the 
currently available data. 

The determination of the Wilson LECs can help to quantify the size of 0(a 2 ) 
terms in a given lattice action. Knowing these LECs, in particular with better 
precision, can significantly contribute to design a lattice fermion action with 
small lattice artifacts, thus allowing to reach the continuum limit in a better 
controlled way. In addition, an independent calculation of the Wilson LECs, 
as carried through here, can be used in chiral perturbation theory fits of light 
meson observables by constraining these fits to lattice data. Thus we think 
that the study of the (connected and full) neutral and charged pion masses of 
this work can be beneficial for many other groups working with Wilson-like 
lattice fermions. 
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